FIRST ORDER DEFORMATIONS OF THE FOURIER MATRIX 

TEODOR BANICA 

Abstract. The N x N complex Hadamard matrices form a real algebraic manifold Cn- 
The singularity at a point H € Cn is described by a filtration of cones T^Cn C T^Cn C 
TjjCn C ThCn, coming from the trivial, affine, smooth and first order deformations. We 
investigate here these cones in the case where H = Fn is the Fourier matrix, i?y = vf 3 
with w = e 2m ' N . Our main result is a simple description of the enveloping tangent space 
TrCn, with the corresponding vectors appearing as plain sums of certain variables. 
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Introduction 

An Hadamard matrix is a square matrix H £ Mjy(±l), whose rows are pairwise or- 
thogonal. The size of such a matrix must be N = 2 or N £ 4N, and the celebrated 
Hadamard Conjecture (HC) states that for any £ 4N, there exists an Hadamard ma- 
trix of order N. Also famous is the Circulant Hadamard Conjecture (CHC), stating that 
the only examples of circulant Hadamard matrices can appear at = 4. We refer to the 
monographs [T5], for a discussion of these key problems. 

There are several interesting generalizations of the Hadamard matrices, see [J]. We will 
be interested here in the complex Hadamard matrices: the square matrices H £ Mjv(C) 
whose entries are on the unit circle, \H%j\ = 1, and whose rows are pairwise orthogonal. 



2000 Mathematics Subject Classification. 05B20 (14B05, 46L37). 
Key words and phrases. Complex Hadamard matrix, Fourier matrix. 

1 



2 



TEODOR BANICA 



The basic example of such a matrix is the Fourier one, P/v = (w^) with w = e 2m ^ N : 

/ 1 1 1 ... 1 \ 

1 w w 2 ... w N ~ l 



N 



In general, the theory of complex Hadamard matrices can be regarded as a "non- 
standard" branch of discrete Fourier analysis. As an example of what can be done with 
the help of such a matrix, we refer for instance to Tao's paper [381. For some other 
potential applications, mainly to quantum physics questions, see [7], [35] , [ID] . 

The complex Hadamard matrices have as well deep connections with von Neumann 
algebras, subfactor theory, knot invariants, spin models, planar algebras, quantum groups 
and free probability theory. See [2], [17], [19], [20], [21], [28], [32]. 

One key problem in the area, raised by Jones in [20], concerns the computation of 
the numbers Ck = dimP^, where P = (P^) is the planar algebra associated to H. These 
numbers, called quantum invariants of H, are in general extremely hard to compute. Very 
little is known about them, and only a few results, conjectures, or even clear questions 
regarding them are available, see [2J, [D], [2D], j3TJ. In particular, we have: 

Problem. What is the relation between the quantum invariants of H , and the geometry 
of the complex Hadamard matrix manifold around H ? 

In order to further comment on this problem, let us recall what the known results on 
the complex Hadamard matrix manifold Cm are. First, this manifold is given by: 

C N = M N (T) n ViVPjv 

Thus, C*tv is an intersection of real algebraic smooth manifolds, so it is a real algebraic 
manifold, not smooth in general. The results about it can be divided into four classes: 

(1) Local structure. The main problem here is of course to describe the singularity of 
Cjv at a given point H e Cn- See pQ, [5], [36], [37] . 

(2) Global structure. At iV < 5 we only have the matrix Pjy and its deformations, as 
shown in [15J. The structure of C% is not known, see [2j, [6], [23] , [34J. 

(3) Rational points. These were introduced by Butson in [TDJ. The main problem here 
is to formulate an analogue of the HC. See [2J, [21], [25j . 

(4) Circulant points. These were introduced by Bjorck in |8j. The main problem here 
is to formulate the correct extension of the CHC. See [3J, [IB] . [27] . 

The present paper mainly belongs to the first direction mentioned above. We will first 
perform some axiomatization work, by introducing the following filtration of cones: 

T^Cn C T^jCn C T h Cn C T h Cn 

Here ThCn is the tangent cone, and T^Cn is the affine tangent cone, obtained by 
restricting attention to the affine deformations. These cones live in the enveloping tangent 
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space T h Cn, obtained by interesecting the tangent spaces to the smooth manifolds M N (T) 
and \fNUNi and contain the trivial tangent cone T^Cn, consisting of vectors which are 
tangent to the trivial deformations, obtained by multiplying rows and columns. 

In general, the computation of these cones is a quite difficult linear algebra problem. 
In the Fourier matrix case, however, we have the following key formula: 

tim(T H C N ) = #{(i,j)\H ij = l} 

This result was established by Tadej and Zyczkowski in [37] , with parts of it going back 
to Karabegov's paper [22], and with the general case discussed in detail in [TJ. 

In this paper, following a number of supplementary ideas from [5], [T3], [35], we will 
obtain a finer result about H = Fn, directly in terms of ThCn'- 

Theorem. For H = Fn the vectors A £ ThCn appear as plain sums of type 

GxHcZ N 

where the L variables are free, and form dephased matrices L GH £ Mgxh(M)- 

We refer to section 5 below for the precise formulation of this result, which requires a 
number of preliminaries, not to be explained in detail right now. 

Let us go back now to the quantum algebra/algebraic geometry problem formulated 
above. This problem was recently investigated in [JJ, with a proposal there involving 
Diaconis-Shahshahani type variables [12]. The various results obtained here, which sub- 
stantially refine those in [JJ, suggest an upgrade of the whole discussion, as follows: 

Question. What algebraic and geometric invariants of H are encoded by the statistics of 
the number of 1 entries, over the equivalence class of H? 

More precisely, consider for instance a usual Hadamard matrix, H £ Mjy(±l). This 
matrix is of course described by the set of indices E C {1, . . . , N} x {1, . . . , N} telling us 
where the 1 entries are. What we propose here is a more geometric approach to H, by 
considering the function (p : Z^ x Z^ — > N which counts the number of 1 entries of the 
various conjugates of H, obtained by switching signs on rows and columns: 

<p(a,b) = #{(i,j)\aibjHij = 1} 

This construction can be generalized to the Butson matrices, and our claim is that (p, 
or just its probabilistic distribution fi, should encode important information about H. 

In general, the computation of fi is a difficult problem. In the real case, this is related 
to the Gale-Berlekamp game [H], [33] and to various questions regarding the 0-1 matrices, 
cf. [IT], [29], [39]. We intend to come back to these questions in some future work. 

The paper is organized as follows: 1-2 are preliminary sections, containing as well a 
number of new results, in 3-4 we state and prove some technical results, and 5-6 contain 
the main result, the probabilistic speculations, and a few concluding remarks. 
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1. Hadamard matrices 

We consider in this paper various square matrices over the complex numbers, H G 
Mjv(C). The indices of our matrices will usually range in the set {0, 1, . . . , N — 1}. 

Definition 1.1. A complex Hadamard matrix is a matrix H G Mn(C) whose entries are 
on the unit circle, \Hij\ = 1, and whose rows are pairwise orthogonal. 

It follows from linear algebra that the columns of H are pairwise orthogonal as well. In 
fact, a matrix H G Mn(T) is Hadamard if and only if U = H/y/N is unitary. 



The basic example is the Fourier matrix F A 



N 



W 



with w = e 2m l N . In matrix 



notation, and by using the above convention for the indices, we have: 





(I 


1 


1 


■ 1 \ 




1 


w 


w 2 


. 


F N = 












V 1 


W N-1 


W 2(N-1) 





There are many other complex Hadamard matrices, see [36]. One simple way of con- 
structing new examples is by taking tensor products, (H ® K) ia j b = HijK ab . 

Observe that the above matrix is nothing but the matrix of the discrete Fourier 
transform, over the cyclic group Z^r. In fact, we have the following result: 



Proposition 1.2. The Fourier matrix of a finite abelian group G — 
the complex Hadamard matrix Fq = F^ 1 <g) . . . <g) F^ k ■ 



x Z 



N k IS 



Proof. For a product of groups G = G' x G" we have Fq = Fc ® Fq", and together with 



the above observation regarding Fjf, this gives the equality in the statement. 

As an example, for the Klein group Z2 x Z2 we obtain the following matrix: 

/ll 1 1\ 
/1 1 \ /1 1 \ i_i 1 

F 2 ® F 2 



□ 



1 



1 

-1 
-1 



-1 
-1 
1 



} -1 

\l -1 -1 I J 
The following notions will play a key role in what follows: 

Definition 1.3. Let H,K G Mn(C) be two complex Hadamard matrices. 

(1) H is called dephased if its first row and column consist of 1 entries only. 

(2) H, K are called equivalent if one can pass from one to the other by permuting rows 
and columns, or by multiplying the rows and columns by numbers in T. 
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In other words, we use the equivalence relation on the N x N matrices coming from 
the action T {A)B) (H) = AHB* of the group G = ((S N x T) x (S N x T))/T. Observe that 
any complex Hadamard matrix can be assumed, up to equivalence, to be dephased. 

We know from Proposition 1.2 that the class of Fourier matrices is closed under tensor 
products. In order to escape from it, we can use deformations. As a basic example, given 
a matrix Q £ M2CT), written Q = (® b d ), we can form the following matrix: 

fa a b b \ 

c —c d —d 



1 

-1 



a b 
c d 



1 

-1 



a a —b—b 
\^c —c —d d J 

With the same data in hand, we can form as well the following matrix: 

fab a b \ 



1 



a b 
c d 



a 

a —b a 
c d —c 
c —d —c 



-b 
-d 
d J 



Observe that the above two matrices are indeed Hadamard. In fact, these matrices 
appear as particular cases of the following general construction, due to Di^a [13] : 

Proposition 1.4. If H £ M/v(C) and K £ Mm{C) are Hadamard, then so are: 

(1) H Q ®K = (QajHijKab) ia,jb, With Q G Mjv/xAf(T). 

(2) H® Q K = (Q lb H tJ K ab ) iaJb , with Q £ M NxM (T). 

These two matrices will be called left and right Dija deformations of H ® K . 

Proof. If we denote by R\ a the ia-th row of H K, we have: 

< R] a , R\ c >= QajHijKa, ■ Q cj H kj K cb = MS ac ^ H l3 H kj = MN5 ac 5 ik 

jb j 

Also, if we denote by Rf a the ia-th row of H ®q K, we have: 

< R 2 ia , R\ c >= QibH %3 K ab ■ Q kb H kj K cb = N5 ik K ab K cb = NM6 ik 5 ac 

jb b 

Thus in both cases we have < Ri a , R kc >= NMd~ ia ^ kc , which gives the result. □ 

Let us take now a closer look at all the above examples of complex Hadamard matrices. 
We can see that the roots of unity often play a key role, and we have: 

Definition 1.5. A complex Hadamard matrix H £ M^(C) is called: 

(1) Of Buston type, if all its entries are roots of unity, of a given order s < 00. 

(2) Regular, if all scalar products between rows decompose as sums of cycles. 
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Here by "cycle" we mean a full sum of roots of unity, possibly rotated by a scalar: 

n 

C = J2 Xe 2kni/n 

k=l 

As a basic example, all the Fourier matrices are of Butson type, and their Di^a de- 
formations are regular. One interest in the above notions comes from the fact that the 
regular complex Hadamard matrices can be fully classified up to iV = 6. See [2]. 

In general, the regularity condition is Definition 1.5 is quite poorly understood: 

Conjecture 1.6 (Regularity). Any Butson matrix is regular. 

More precisely, for exponents of type s = p a or s = p a q b one can prove that any 
vanishing sum of s-roots of unity decomposes a sum of cycles, but this is not true in 
general. See [23]. Here is actually the simplest counterexample, with w = e 2? ™/ 30 : 

w 5 + w 6 + w 12 + w 18 + w 24 + w 25 = 

So, what Conjecture 1.6 says is that such a "tricky sum" cannot be used for constructing 
a complex Hadamard matrix. We refer to [2] for further details on this conjecture. 

Let Cjv be the real algebraic manifold formed by the NxN complex Hadamard matrices. 
We denote by M x an unspecified neighborhood of a point in a manifold, x G M. Also, for 
q G Ti, meaning that q G T is close to 1, we define q r with r G M by (e lt ) r = e ltr . 

Proposition 1.7. For H G Cjv and A G Mn(M>), the following are equivalent: 

(1) = Hijq Aij is an Hadamard matrix, for any g 6 Ti. 

(2) Ylk HikHjkq Aik ~ Ajk = 0, for any % ^ j and any q G Ti. 

( 3 ) J2k Hi k H jk ip(A ik - A jk ) = 0, for any i ^ j and any (/):!-)•£. 

( 4 ) J2kGE r H ik H jk = for any i ^ j and r G R, where E r M = {k\A ik - A jk = r}. 

Proof. All the equivalences are elementary, and can be proved as follows: 

(1) <^=^ (2) Indeed, the scalar products between the rows of H q are: 

< H?,H] >= J2 H ikq Alk H Jk q Ajk = ^#^ fc g Alfc ^ fc 

k k 

(2) (4) This follows from the following formula, and from the fact that the power 
functions {q r \r G M} over the unit circle T are linearly independent: 

Y,H lk H jk q A ^- A ^ = J2v r Yl H ^ k 

k reM keET. 

(4) =^> (3) This follows from the following formula: 

^2H ik H jk Lp(A ik - A jk ) = s ^ip(r) 22 H ik H jk 

k reM k&ET. 

(3) (2) This simply follows by taking ip(r) = q r . □ 
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Observe that in the above statement the condition (4) is purely combinatorial. 

In order to understand the above type of deformations, it is convenient to enlarge 
attention to all types of deformations. We keep using the neighborhood notation M x 
introduced above, and we consider functions of type / : M x — > N y , which by definition 
satisfy f(x) = y. With these conventions, we introduce the following notions: 

Definition 1.8. Let H G Mw(C) be a complex Hadamard matrix. 

(1) A deformation of H is a smooth function f : Ti — > (Cn)h- 

(2) The deformation is called "affine" if fij(q) = Rijq Aii , with A G M N (M). 

(3) We call "trivial" the deformations of type fij{q) = Hijq ai+b] , with a, b G M. N . 



Here the adjective "affine" comes from /« (e l 



because the function t — > A^t 
which produces the exponent is indeed affine. As for the adjective "trivial", this comes 
from the fact that f(q) = (H i jq ai+bj ) i j is obtained from H by multiplying the rows and 
columns by certain numbers in T, so it is automatically Hadamard. See [36J. 

The basic example of an affine deformation comes from the Di^a deformations H ®q K, 
by taking all parameters G T to be powers of q G T. As an example, here are the 
exponent matrices coming from the left and right Dita deformations of F2 (g> F2. 

(a a b b\ ( a b a b\ 



A, 



d 
b 
d 



d 
b 

dj 



b 
d 
d 



b 
d 

dj 



In order to investigate the above types of deformations, we will of course use the cor- 
responding tangent vectors. So, let us first recall that the manifold Cat is given by: 

C N = Mjv(T) n Vnu n 

This observation leads to the following definitions, where in the first part we denote by 
T X M the tangent space to a point in a smooth manifold, x G M: 

Definition 1.9. Associated to a point H G CV are the following objects: 

(1) The enveloping tangent space: T h Cn = T H M N (T) fl T h \/NUn. 

(2) The tangent cone ThCn: the set of tangent vectors to the deformations of H. 

(3) The affine tangent cone T^C^: same as above, using affine deformations only. 

(4) The trivial tangent cone T^Cn: as above, using trivial deformations only. 

Observe that ThCn,T^Cn are real linear spaces, and that ThCn,T^Cn are two-sided 
cones, in the sense that they satisfy A G M, A G T XA G T. 
Observe also that we have inclusions of cones, as follows: 

T^Cn C T^jCm C T h Cn C T h Cn 

In more algebraic terms now, the above tangent cones are best described by the corre- 
sponding matrices, as follows: 
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Theorem 1.10. The cones T^C^ C T^C^ C T H C^ C T H C^ are as follows: 

(1) ThCn can be identified with the linear space formed by the matrices A G Mw(R) 
satisfying J2 k H ik H jk (A ik - A jk ) = 0, for any i,j. 

(2) ThCn consists of those matrices A £ Mjy(R) appearing as A^ = g^(0) ; where 
g : M N {R) ->■ Mjv(R) safe/^es J2k H ik H jke i{9ik{t) ~ 9jk{t)) = /or any z,j. 

(3) T^Cjv is formed by the matrices A G Mjv(R) satisfying s ^ jk H ik Hj k q Aik ~ A i k = 0, 
for any i ^ j and any q G T. 

(4) T^Cn is formed by the matrices A £ Mjy(R) which are of the form A^ = a; + 6j ; 
/or certain vectors a, b G R . 

Proof. All these assertions can be deduced by using basic differential geometry: 

(1) This result is from [T], the idea being as follows. First, M^(T) is defined by the 
algebraic relations \Hij\ 2 = 1, and with = JQj + iY^ we have: 

'I 'I, 2 = '/S AO + ^) = 2(X«X« + Y.Y.j) 

Now since an arbitrary vector £ G ThM^(C), written as £ = . acijXij+PijYij, belongs 
to TffMjv(T) if and only if < £, d|ii/y| 2 >= for any i,j, we obtain: 



t h m n (t) = { r^^i^-x^ 




G 1 

We also know that VNUn is defined by the algebraic relations < Hi,Hj >= NSij, 
where Hi, . . . , Hn are the rows of H. The relations < Hi, Hi >= N being automatic for 
the matrices H G Mjv(T), if for % ^ j we let =< i^j, iifj >, then we have: 

T^CV = {£ G T H M N (T)\ < £, 4 >= 0, Vi ^ i} 
On the other hand, differentiating the formula of L^- gives: 

^ij = ^^(Ajfc + — ii^-jfc) + (Xjk — iYj k )(X ik + iY ik ) 

k 

Now if we pick £ G TjjM^(T), written as above in terms of A G Mjy(R), we obtain: 

< C, >= « H ik Hj k (A ik — Aj k ) 

k 

Thus we have reached to the description of ThCn in the statement. 
(2) Pick an arbitrary deformation, and write it as fij(e lt ) = H^e 19 **®. Observe first 
that the Hadamard condition corresponds to the equations in the statement, namely: 

J2Hi k H jk e l ^ kit) -^ k{t)) = 
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Observe also that by differentiating this formula at t = 0, we obtain: 

k 

Thus the matrix Ay = 5^(0) belongs indeed to ThCn, so we obtain in this way a 

certain map ThCn — > ThCn- In order to check that this map is indeed the correct one, 
we have to verify that, for any i, j, the tangent vector to our deformation is given by: 

Q; — 9ij(0)(YijXij — XijYij) 

But this latter verification is just a one- variable problem. So, by dropping all i, j indices 
(which is the same as assuming N = 1), we have to check that for any point H e T, written 
H = X + iY, the tangent vector to the deformation f(e lt ) = He 19 ^ is: 

£ = g\0)(YX-XY) 

But this is clear, because the unit tangent vector at H e T is rj = — i(YX — XY), and 
its coefficient coming from the deformation is (e* 9( *))j t=0 = — ig'(0). 

(3) Observe first that by taking the derivative at q — 1 of the condition (2) in Propo- 
sition 1.7, of just by using the condition (3) there with the function <p(r) = r, we get: 

H ik Hj k (p(A ik — A jk ) = 

k 

Thus we have a map T^Cn — > ThCn, and the fact that is map is indeed the correct 
one comes for instance from the computation in (2), with gij(t) = Aijt. 

(4) Observe first that the Hadamard matrix condition is satisfied: 

J2 HikHj k q Aik ~ Ajk = q a ^ ^ H ik H jk = 5 tJ 

k k 

As for the fact that T h Cn is indeed the space in the statement, this is clear. □ 

Let now Dn C Cn be the real algebraic manifold formed by all the dephased N x N 
complex Hadamard matrices. Observe that we have a quotient map Cn — > Dn, obtained 
by dephasing. With this notation, we have the following refinement of (4) above: 

Proposition 1.11. We have a direct sum decomposition of cones 

T h Cn = T h Cn © T^Dn 

where at right we have the affine tangent cone to the dephased manifold Cn — > Dn- 

Proof. As explained in [TJ, if we denote by M N (M) the set of matrices having outside 
the first row and column, we have a direct sum decomposition, as follows: 

T H C N = M N (R)®T H D N 
Now by looking at the affine cones, and using Theorem 1.10, this gives the result. □ 
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2. The defect, revisited 

The following key definition, whose origins go back to the work of Karabegov [22] and 
Nicoara [30], was given by Tadej and Zyczkowski in [37J: 

Definition 2.1. The defect of a complex Hadamard matrix H e Cn is the real dimension 
d(H) of the enveloping tangent space ThCn = ThMn(T) fl ThVNUn- 

We should mention that there are actually several notions of defect, and terminologies 
for them, see [5], [37]. In what follows we will use the above notation and terminology. 
Here few basic properties of the defect: 

Proposition 2.2. Let H G Cn be a complex Hadamard matrix. 

(1) IfH~H then d(H) = d(H). 

(2) We have 2N — 1 < d(H) < N 2 . 

(3) If d(H) = 2N — 1, the image of H in the dephased manifold Cn — >■ F>n is isolated. 

Proof. All these results are well-known, the proof being as follows: 

(1) Consider indeed the equations for the enveloping tangent space, namely: 

^^H ik Hj k (A ik — Ajk) = 
k 

If we let now = aibjHij with \a,i\ = = 1 be a trivial deformation of our matrix 
H, the equations for the enveloping tangent space for K are: 

^^aibkHikdjbkHjk(Aik — Ajk) = 
k 

By simplifying we obtain the equations for H, so d(H) is invariant under trivial defor- 
mations. Since d(H) is invariant as well by permuting rows or columns, we are done. 

(2) Consider the inclusions T^Cn C ThCn C ThCn- Since dim(T^CV) = 2iV — 1, the 
inequality at left holds indeed. As for the inequality at right, this is clear. 

(3) In case of equality in the inequality at left in (2), we must have ThCn = T^Cn- 
Thus in case of equality any deformation of H is trivial, and it follows that the image of 
H in the quotient manifold Cn ~ > Dn is indeed isolated, as stated. □ 

As a first remark, the upper bound in Proposition 2.2 (2) is of course trivial. We believe 
that the inequality d{H) < jV( - j ^ +1 ^ should hold in general, but we don't have a proof for 
this fact. We will be back to this issue a bit later, after Theorem 2.6 below. 

In the case of Fourier matrices, the computation of the defect is as follows: 

Proposition 2.3. Let F be the Fourier matrix of a group G = Zjvi x . . . %N k - 

(1) TfCn = \PF*\Pij = Pi+j.j = Pi-j}- 

(2) d(F) = J2 gEG [G:<g>]. 

(3) d(F) is also the number of 1 entries of F. 
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Proof. This is basically known from [22], [37], with some slight improvements and gener- 
alizations coming from [1], [5], and the proof is as follows: 

(1) According to Theorem 1.10 (1), the equations are J2k FikFjk{Aik — Ajk) = 0. But 
these equations simply read (AF)^^ — (AF)j ti ^j = 0, and together with the fact that 
A = (AF)F* must be real, this gives the conditions in the statement. See [TJ. 

(2) The proof here uses an identification of real vector spaces, as follows: 

{P e M N (C)\P i:j = P i+jjj = iVJ < g >,R) 

g&G 

Indeed, if we let G 2 = {g G G\2g = 0}, and then choose a partition of type G = 
GjLllU (— X), the formula giving the above identification is P = ®P g , with: 

f Pa (J e G 2 ) 

P g (t) = I Re(P^ (j G X) 

\lmiPy) (je-X) 
With this identification in hand, the result follows from (1). 

(3) This observation, from [22J, follows from (2), and from the definition of F. □ 

It is perhaps not obvious from Proposition 2.3 that the matrices A = PF* constructed 
in (1) there are real. Here is an equivalent formulation of the result: 

Proposition 2.4. For F = Fa, the space TpC^ consists of matrices of the form 

A ik = Re 



J2 wfk Q^ + Im \J2 wfkR ^ 



where Q, R G Mg(R) satisfy Qij = Qi+jj = Qi-j and Rtj = Ri+jj = —Ri 



Proof. Indeed, let us pick P G Mg(C), and write it as P = Q + iR, with Q,Re Mq{ 
Then the conditions Pjj = Pi+jj = P-j in Proposition 2.3 (1) are equivalent to the 
conditions on Q, R in the statement, and the matrix A = PF* is given by: 

A ik = iJ2 w ' fk Qif + J2 w ~ fkRi f 

k f 

= -\ 5> /fc - w ~ fk )Q« + \ 5> /fc + w ' fk )^f 
f f 

But this gives the formula in the statement, and we are done. □ 

Proposition 2.5. The defect of basic matrices of type F = Fq is as follows: 

(1) For G = Z N with N = [7;Pf we have d(F) = N]]^! + a { , - |). 

(2) In the Walsh matrix case G = 7Lt 1 we have d[F) = N % , with N = 2 k . 
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Proof. We use Proposition 2.3, and more specifically the assertion (3) there: 

(1) This formula follows indeed from Proposition 2.3. See [37] . 

(2) Observe first that the matrix in the statement is indeed the k-th Walsh matrix, 
W 2k = F® k . We recall now that, according to our tensor product notations: 

fW 2k -i W 2 k-i \ 
2k \W 2 k-i -W 2 k-iJ 

Now if we denote by 4 the defect of W 2 k , which is the number of 1 entries, then: 

4 = 34-i + (4*- 1 - 4-i) = 24-i + 4 fe - 1 

By solving the recurrence we obtain 4 — 2 fe ~ 1 (2 fc + 1), which gives the result. □ 

Let us try now to better understand the formula d(F) = jV ^ +1 - > found in Proposition 
2.5 (2) above. For this purpose, best is to get back to Theorem 1.10 and its proof. 
Everything there was based on the formula CV = M N (T) fl VNUn, and more specifically, 
on the idea of performing all computations inside the space ThM^{T). But, we can of 
course try to perform as well all computations inside the space ThVNUn. 

Let us denote by Mn(C) sci the iV 2 -dimensional real vector space formed by the complex 
N x N matrices which are self-adjoint, in the sense that they satisfy M = M*. 

Theorem 2.6. Let H G Cn be a complex Hadamard matrix. 

(1) We have f H C N ~ {S G M JV (C) sa |(5iJ) ii i? ij G K,Vz,j}. 

(2) The correspondence A — )■ S is given by S = MH* , with = A^H^. 

(3) For H = F G we have f F C N ~ {S G M^CHS^ = S t+jA }. 

(4) For Walsh matrices, H = W N with N = 2 k , we have f H C N = M N (R) symm . 

Proof. The idea is that (1,2) follow by carefully rewriting the proof of Theorem 1.10, and 
that (3,4) are consequences of (1,2). More precisely: 

(1) Since \/NUn is defined by Ly = N5ij, with Ly =< Hi, Hj >, we have: 

T H ^NU N = {£ G T H M N (C)\ < >= 0,Vz, j} 
On the other hand, the derivative of Lij =< Hi,Hj > is given by: 

Lij = ^^(X ifc + iY ik )(Xj k — iYjk) + {Xjk — iYj k )(X ik + iY ik ) 

k 

Now with £ G T H M N (C), written as £ = + PiYiji we have: 

< £, L i:j >= ^(«jfc + i(3 jk )(X ik - iY ik ) + (a ik - i(3 ik )(X jk + iY jk ) 

k 

Thus, in terms of the complex numbers jij = aij + i/3y, we obtain: 
< £, Ly >= ^2-f jk H ik + ji k H jk = (jH*)ji + {^H*)ij 

k 
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We therefore conclude that for £ = £\ a-iX^ + faY^ as above, we have: 

£ G T H M N (C) 1 H* G M N (C) sh 

Here M]y(C) sh is the set of N x N complex matrices which are skew-hermitian, in the 
sense that they satisfy M* = —M. Now recall from the proof of Theorem 1.10 that we 
have £ G T H M N (T) if and only if the following scalar products vanish, for any 

< £, XijXij + YijYij >= <\;,.Xjj + PijYij = Re((aij + ip ij )(X ij - iY i:j )) = lit (' 

Thus we have reached to the following description of the enveloping tangent space: 

f H C N — {7 G M^C^H* G M N (C) sh , ll:i H tJ G iR,Vi, 3} 

In terms of the matrices S = i'yH*, we therefore obtain, as claimed: 

f H C N — {S G M^CrKSH)^ G R,Vi, j} 

(2) Let us pick a matrix A G Mjv(R) as in Theorem 1.10 (1), and consider the associated 
vector £ G T H C N , given by £ = £V A^^X^ - X^Y^). We have then: 

Im^'jijHij) = ftijXij — (XijYij = —AijX^j — A^Y^ = —Aij 

Thus the bijection A — )• 7 we are looking for is given by 7^-^^ = —iA^, which reads 
7jj = —iAijHij. In terms of the above matrices S* = i'jH*, this shows that the corre- 
sponding bijection A — > S is given by S = (AijHij)ijH* , and we are done. 

(3) Let P G Mjv(C) be as in Proposition 2.3. Then A = PF* is given by: 

^ = (PF*) - = J>-*''P ifc 

k 

Let us compute now the matrices M, 5 in (2) above. First, we have: 

M ij = A ij F ij = Y t ™ ( ' i ~ k)j Pik 
k 

Thus the matrix S = MF* is given by: 

S u = (MF*) U = J2M lJ (F*), l = Y,^ l ~ k ~ l)] P lk = Pi,i-i 

j jk 

So, we have P i3 = and the conditions Pij = Pi+jj = Pi-j become: 

Q. . . — Q. . . — Q. . . 

U i,i-J ^1,1+3 

Now observe that the equality on the right simply tells us that S is self-adjoint. Thus, 
what is left is the equality on the left, which gives the formula in the statement. 

(4) Since Wn is real, so are the corresponding matrices S, and the result follows from 
(1), because the condition there on S is automatically satisfied. □ 
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Observe that Theorem 2.6 provides a more conceptual proof for the defect formula 
d(Wjsr) = N ( N 2 +1 ^ in Proposition 2.5. As already mentioned, we believe that the inequality 

d{H) < jV ^ +1 - > should hold in general, but we don't have a proof for this fact. 
Let us discuss now some rationality questions: 

Definition 2.7. The rational defect of H £ Cn is the following number: 

d Q (H) = dim Q (T H C N n Mjv(Q)) 

The vector space on the right will be called rational enveloping tangent space at H. 

As a first observation, this notion can be extended to all the tangent cones at H, and 
by using an arbitrary field KcC instead of Q. Indeed, we can set: 

T* H C N (K)=T* H C N nM N (K) 

However, in what follows we will be interested only in the objects constructed in Defi- 
nition 2.7. It follows from definitions that d<q(H) < d(H), and we have: 

Conjecture 2.8 (Rationality). For Butson matrices we have dq(H) = d(H). 

In fact, the original statement in [I] is that the above equality should hold in the regular 
case. However, since the regular case is not known to fully cover the Butson matrix case, 
as explained in Conjecture 1.6, we prefer to state our conjecture as above. 

Let Ctv(s) be the set of N x N complex Hadamard matrices having as entries the s-th 
roots of unity. With this notation, we have the following elementary result: 

Proposition 2.9. The rationality conjecture holds for H £ Cn{s) with s = 2, 3, 4, 6. 

Proof. Let us recall that the equations for the enveloping tangent space are: 

/ y HikHjk(A ik — A jk ) = 

k 

In the case s = 2 these equations are all real, and have rational (±1) coefficients. In 
the case s = 3, 6 we can use the well-known fact that, with j = e 2m ^ 3 , the real solutions 
of x + jy + j 2 z = are those satisfying x = y = z, and we conclude that once again 
our system, after some manipulations, is equivalent to a real system having rational 
coefficients. As for the case s = 4, here the coefficients are —1, —i, so by taking the 
real and imaginary parts, we reach once again to system with rational coefficients. 

Thus, in all cases under investigation, s = 2, 3, 4, 6, we have a real system with rational 
coefficients, and the result follows from standard linear algebra. □ 

Observe that the above method cannot work at s = 5, where the equation a + wb + 
w 2 c + w 3 d + u> 4 e = with w = e 27 ™/ 5 and a, b, c, d, e £ K can have "exotic" solutions. 

We will prove in section 5 that Conjecture 2.8 is verified for the Fourier matrices. 

Finally, let us go back to Proposition 2.3, and to the formula d(H) — |1 £ H\ found 
there. We expect this formula to be valid under much more general assumptions: 
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Problem 2.10. When does the defect formula d(H) — |1 6 H\ hold? 

We believe for instance that this formula should hold for matrices of type H = Fg®qFk, 
under fairly general assumptions on the deformation matrix Q 6 Mjvxm(T). For instance 
Proposition 3.5 below suggests that the result should hold for Q dephased and elsewhere 
generic. However, we do not know what the exact assumptions on Q should be. 

3. Tensor product gluing 

In this section we discuss the computation of the various tangent cones at a tensor 
product H £g> K. The problem that we are interested in is to determine how the various 
tangent cones at H, K, from Definition 1.9 above, glue together at H (g) K. 

We recall that our tensor product convention is (H (g) K) ia j b = HijK ab . In matrix form, 
by using the lexicographic order on the double indices: 



H U K ... H 1N K 



H N1 K ... H NN K, 
Let us begin with a elementary result, regarding the enveloping tangent spaces: 

Proposition 3.1. We have T H C N ® T K C M C T H ® K C NM - 

Proof. Indeed, for a matrix A = B <8> C, we have the following formula: 

ia,kc 

(H ® K) jb kc A iatkc = 22 HikK ac Hj k K bc B ik C ac 



kc kc 



Hi k Hj k B ik K ac K bc C a 



On the other hand, we have as well the following formula: 



^2(H ®K)ia >kc (H ® K) jbM A jb)kc = ^2H ik K ac H jk K bc B jk C bc 

kc kc 

= H ik Hj k Bj k K ac K bc C bc 

k c 

Now by assuming B e T H C^ and C G T k Cm, the two quantities appearing above on 
the right are equal. Thus we have indeed A e Th®kCnm, and we are done. □ 

In general, it is not clear when the inclusion in Proposition 3.1 is an isomorphism. We 
refer to |35J and to section 5 below for a number of results in this direction. 

In order to further study the tangent cones, let us first review in detail the notion of 
Dita deformation, from Proposition 1.4 above. We recall from there that: 

H q® K = (QajHijKab)ia,jb, H ® Q K = (Qi b HijK ab ) ia:jb 
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In matrix notation, by using more convenient lower-case letters for the entries of Q, 
and by using the lexicographic order on the double indices, we have: 

QmH 1N Kn . . . q 1N H 1N KiM \ 



H, 



K = 



( quHuKu 
QmiHhK M i 

QuHniKh 
\QmiHn\Kmi 



qiiHnKiM 
QmiHhK mm 



QmnHinK M i . . . q MN Hi N K MM 



QmnHnnKmi 



QmnHnnKmm 



QmiHniKmm 

As for the right Dfja deformation, the entries here are of course the same as above, but 
up to a certain permutation. More precisely, the formula is: 



H® Q K 



( quHuKu . 

QniHniKu . 
\c[niHniKmi ■ 



QimHuKim 
QimHuK mm 



QnmH N iK 



1M 



QniHnnKh . 



QimHinKi M \ 
QimHinK mm 



QnmHnnKim 
InmHnnKmm 



QnmHniKmm 

The fact that the entries of the right and left Di^a deformations are the same modulo 
a permutation is in fact a consequence of the following simple fact: 

Proposition 3.2. We have an equivalence H q® K ~ K ®q H. 

Proof. According to the formulae in Proposition 1.4 above, we have: 

(H q® K) ia j b = q aj HijK ab = q aj K ab Hij = (K <g)<g H) aijbj 

Now since the transformation M ia j b — > M ai ^ b j is implemented by certain permutations 
of the rows and columns, the above two matrices are indeed equivalent. □ 

Observe now that, if we look at the complex Hadamard matrices modulo the equiv- 
alence relation in Definition 1.3, we can always assume that our parameter matrix Q is 
"dephased" , in the sense that its first row and column consist of 1 entries only. 

As an illustration here, at iV = M = 2 we have the following result: 



Proposition 3.3. Any Di$a deformation o/F 2i2 



2,2 



1 1 

1 q 



F 2 <8> F 2 is equivalent to 

/l 1 1 l\ 

-1 1 -1 

q -1 -q 

-q -1 q ) 



1 
1 

V 
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for a certain value of the parameter q G T. 

Proof. First, by using Proposition 3.2 we may restrict attention to the case of right Dita 
deformations. But here, with Q = (% £), our claim is that we have F 2 ®q F 2 ~ F 2 9 2 , with 
q = ad /be. Indeed, by dephasing our matrix we obtain: 



/ a 


b 


a b \ 


(I 


1 


1 1 \ 




(I 


1 1 


i \ 


a 


-b 


a —b 


1 


-1 


1 -1 




1 


-1 1 


-l 


c 


d 


— c —d 


- c 
a 


d 


c d 
a b 




1 


ad i 

ad i 
be 


ad 


V 


-d 


—c d J 


\ a 


~ b d 

b 


c d 
a b / 




V 


ad 

be / 



But this means F 2 ®q F 2 — F 22 i with q = ad/bc, and this finishes the proof. □ 

Let us discuss now the computation of the various tangent cones for a Dita deformation. 
Here we basically have just one result, concerning the case where the deformation matrix 
is dephased and elsewhere generic, in the following sense: 

Definition 3.4. A rectangular matrix Q G Mn x mO^) is called "dephased and elsewhere 
generic" if the entries on its first row and column are all equal to 1, and the remaining 
(N — 1)(M — 1) entries are algebrically independent over Q. 

Here the last condition takes of course into account the fact that the entries of Q 
themselves have modulus 1, the independence assumption being modulo this fact. 
We have the following result, extending the 4x4 computations in [Tj: 

Proposition 3.5. If H G C^,K G Cm are dephased, of Butson type, and Q G M NxM (T) 
is dephased and elsewhere generic, then A = (A ia>kc ) belongs to T h ® qK Cnm iff 

4L = *&, A& = A&, {A%) xy ef K C M 
hold for any a, b, c and i ^ j , where A l J c = J2 k HikHjkA ia ^ c - 
Proof. Consider the system for the enveloping tangent space, namely: 
^2(H ® K) iaM (H ® K) jh kc {A iaM - A jbM ) = 

fee 

We have (H Cg> K) ia j b = q ib HijK ab , and so our system is: 

^ QicQjcK a cK bc ^ HikHjk(Ai aj kc ~ Ajh^kc) = 
c k 

Consider now the variables A l j c = ^2 k HikHjkA iat k c in the statement. We have: 

Aac = ^ HikHjkAiahc = ^ HjkHikAiakc 

k k 
Thus, in terms of these variables, our system becomes simply: 

J2QicqjcK ac K bc (A i i-Aj c ) = 
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More precisely, the above equations must hold for any i,j,a,b. By distinguishing now 
two cases, depending on whether i,j are equal or not, the situation is as follows: 

(1) Case i ^ j. In this case, let us look at the row vector of parameters, namely: 

{qic<ijc)c = (1, QnQji, ■ ■ ■ , QiN<ijN) 

Now since Q was assumed to be dephased and elsewhere generic, and because of our 
assumption i ^ j, the entries of the above vector are linearly independent over Q. But, 
since by linear algebra we can restrict attention to the computation of the solutions over 

Q, the i ^ j part of our system simply becomes A l J c = A 3 b l c , for any a, b, c and any i ^ j. 
Now by making now a, b, c vary, we are led to the following equations: 

Aac = ^-bci ^ac = ^"c, Vd, 6, C, % ^ j 

(2) Case % = j. In this case the parameters cancel, and our equations become: 

K ac K bc {A™ c - Ag) = 0, Va, 6, c, i 

c 

On the other hand, we have A™ c = J2 k A iajkc , and so our equations become: 

Y,K ac K bc {At-Al) = ^ Va,b,c,i 

c 

But these are precisely the equations for the space T k Cm, and we are done. □ 

Let us go back now to the usual tensor product situation, and look at the affine cones. 
In general, we cannot expect T^^ k C^m to appear from a simple gluing procedure, so the 
question which makes sense is rather that of finding the biggest subcone of T^j^ k CnM) 
obtained by gluing TfjC n ,T^C m ■ Our answer here, which takes into account the two 
"semi-trivial" cones coming from the left and right Di^a deformations, is as follows: 

Theorem 3.6. The cones T^Cn = {B} and T^Cm = {C} glue via the formulae 

Aiajb = A-Bjj + IpjCab + X ia + Yj b + F a j 
Aiajb + f^Cab + X ia + Yjb + Eib 

producing in this way two subcones of the affine cone T^^ k Cnm = {A}. 

Proof. Indeed, the idea is that X ia , Yjb are the trivial parameters, and that E ib , F a j are 
the Di^a parameters. In order to prove the result, we use the criterion in Theorem 1.10 
(3) above. So, given a matrix A = (A ia jb), consider the following quantity: 

P = J2 H ikH jk K ac K bc q A ^- A ^ 

kc 
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Let us prove now the first statement, namely that for any choice of matrices B £ 
T^Cn,C G TjjCm and of parameters \,ipj, X ia ,Yj , F a j, the first matrix A = (A ia j ) 
constructed in the statement belongs indeed to T^^ k Cnm- We have: 

A-Bjfc + i>kC ac + X ia + Y kc + F afe 
Now by substracting, we obtain: 

Aia,kc — ^jb,kc = ^{Bik — Bj k ) + i)k{C ac — Cbc) + (X ia — Xj b ) + (F ak — F bk ) 

It follows that the above quantity P is given by: 

p — H uH uK Kk a^ Bik ' B ^ + ^ k{Cac ^ Cb ^ + ^ Xia ^ x ^ +i ~ Fak ' Fbk ^ 

kc 

= q x ia -X jb H ik H jk q F « k - Fbk q x ^ k - B ^ ^ K ac K bc (q i ' k f ac ~ Cbc 

k c 

= 5 ab q Xi "- x i* J2 HikH jk (q X ) Bik - Bik = <U<% 

k 

Thus Theorem 1.10 (3) applies and tells us that we have A e T^^ k Cnm, as claimed. 
In the second case now, the proof is similar. First, we have: 

^-ia,kc = 4>cBik + flC ac + Xi a + Y kc + Efc 

Aj b) kc = 4> c Bjk + ^C bc + Xj b + Y kc + Ej C 
Thus by substracting, we obtain: 

Aia,kc — Aj bjkc = (j) c (B ik — Bj k ) + fi(C ac — C bc ) + (Xia ~ Xj b ) + (E{c — Ejc) 

It follows that the above quantity P is given by: 

p — 'y y H k H k K K q ( t' c ( Bik ~ B 3 k ') + f J, (c : ac-c bc )+(Xi a -x jb )+(Ei C -Ej C ) 

kc 

= q x ™- x *> K ac Kbcq Elc ~ E3C q^ Cac ~ Chc) ^ H ik H jk (q^ c ) B ^ B ^ k 

c k 

= 6 ijq x *- x » K ac K bc ( q n Cac ~ Cbc = SijSai, 

c 

Thus Theorem 1.10 (3) applies again, and gives the result. □ 

We believe Theorem 3.6 above to be "optimal", in the sense that nothing more can 
be said about the affine tangent space T^ k Cnm, in the general case. This claim is 
supported by various computations for F N <g> F M , and by results in [5], [37] . 

In view of these observations, the following problem appears: 

Problem 3.7. What is the affine cone gluing at an iterated tensor product of complex 
Hadamard matrices, Hi <g> . . . <g> H n ? 
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In principle the answer can be simply obtained here by an iterated application of The- 
orem 3.6. However, it is not clear how to best formulate the final result. 

4. Fourier matrices 

In this section and in the next one we study the Fourier matrix F^, and we state and 
prove our main result, concerning the enveloping tangent space at Fn. 

There is a long story here, going back to Karabegov's paper [22], then to the work of 
Dita [13], Tadej-Zyczkowski j3S], [37], and Barros e Sa and Bengtsson [5J. Our result here 
improves the previously known facts about ThCn, and can be viewed as a complement, 
or new ingredient, to the current exploration of the tangent cone ThCn in [5]. 

Since the whole subject belongs to a new and somewhat mysterious area of linear 
algebra, we will give many examples and details. Let us first discuss in detail the isotypic 
case, where N = p a is a prime power, following [5], [36]. First, we have: 

Proposition 4.1. For F = F p with p prime, we have 

TiXp rpo f~1 rp /~t rp r~i 

pLyp — 1 pLy p — 1 H^p ~ 1 F^p 

and this space consists of the matrices of type Aij = aj + bj, with a, 6 6 R p . 

Proof. It is enough to check that we have TpC p = TpC v , which means that any matrix 
A G TpCp should decompose as = + bj, with a, b G MP. At p = 2 we have: 

"=(; '*=0-i 

Thus the exponent matrices follows, with x + t = y + z: 

A = PF * = fa + c a-c\ = fx y 
\b + c b — cj \zt 

This gives the result. At p = 3 now, with w = e 2m ^ 3 we have: 

fa z z\ /ill 
P = lb z z , F* = 1 w 2 
\czzl \l w 

Thus the exponent matrices have a similar look, and this gives the result: 

(a + z + z a + w 2 z + wz a + wz + w 2 z^ 
b + z + z b + w 2 z + wz b + wz + w 2 z 
c + z + z c + w 2 z + wz c + wz + w 2 z 

In general now, we have Aqq + A^ = A^ + Aqj, and this gives the result. □ 

Observe that Proposition 4.1 is somehow not fully satisfactory, because the fact that 
the defect is d(F) = 2p — 1 is not entirely obvious from it. So, let us improve it: 
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Proposition 4.2. For F = F p , the elements A e TpC v are the solutions of 

'01 i rlO 



A - — r 00 



i rui , rll 

00 + + Ho 



where the L variables are free, with Lq] = for j ^ 0, and L}q = for i ^ 0. 
Proof. The claim is that there exist free variables L, such that: 

^oo = 



r00 
-^00 



A 0j = 

A i0 = 
Aij = 



(<^o) 



rOO , r 10 
-^00 -^iO 



rOO i rui , r 



01 



v 00 



10 
v i0 



But this follows from Proposition 4.1, by taking Lqq = A o, and so on. 



□ 



Let us discuss now the case N = p 2 . We recall that M°(R) is by definition the set of 
real n x n matrices having on the first row and column. We will need the following key 
result, regarding the enveloping tangent space to the dephased manifold CV — > D^: 



Proposition 4.3. For F = F p 2 with p prime, we have 

f F D N = {Ae M N {M)\3M e M;(R), Aij = 
where i,j are the reminders ofi,j modulo p. 
Proof. Let us first work out the case p = 2. Here the defect is 



A i i 

1 -i -1 

1 -1 1 

\1 i -1 



I, and: 
1\ 



i 

-1 
—i 



( a z e z\ 

b z f z 

c z e z 

\d z f z J \1 i -1 —i J 

Thus with z = p + iq we obtain the following formula: 

/ a + e + 2p a — e + 2q a + e — 2p a — e — 2g\ 

b + f + 2p b-f + 2q b + f-2p b-f-2q 

c + e + 2p c — e + 2q c + e — 2p c — e — 2q 

\d + f + 2p d-f + 2q d + f-2p d - f - 2q J 

Now by assuming that A is dephased, as in the statement, we obtain that we have 
a = c = e= p = q = and b = d = —f, and so our matrix is given by: 



A = PF* 



A = 



/0 \ 

-2f -2f 



\0 -2f -2fJ 
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Thus we reach to the conclusion in the statement. At p = 3 now, we have: 



fa 


X 


y 


u 


z 


z 


u 


y 


x\ 




n 


1 


1 


1 


1 


1 


1 


1 


1\ 


b 


X 


y 


V 


z 


z 


V 


y 


X 




1 


w 8 


7 

w 


w 6 


w 5 


4 

W 


w 3 


w 2 


w 


c 


X 


y 


w 


z 


z 


w 


y 


X 




1 


7 

w 


w 5 


w 3 


w 


w 8 


w 6 


4 

w 


w 2 


d 


X 


y 


u 


z 


z 


u 


y 


X 




1 


w 6 


w 3 


1 


w 6 


w 3 


1 


w 6 


3 

w 


e 


X 


y 


V 


z 


z 


V 


y 


X 




1 


w 5 


w 


w 6 


w 2 


7 

w 


w 3 


w 8 


4 

w 


f 


X 


y 


w 


z 


z 


w 


y 


X 




1 


4 

w 


w 8 


w 3 


7 

w 


w 2 


w 6 


w 


w 5 


9 


X 


y 


u 


z 


z 


u 


y 


X 




1 


w 3 


w 6 


1 


w 3 


w 6 


1 


w 3 


w 6 


h 


X 


y 


V 


z 


z 


V 


y 


X 




1 


w 2 


4 

w 


w 6 


w 8 
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By assuming now that A is dephased, we obtain the following equations: 

a = x = y = z = u = 0, d = g = 



b = e = h = —(v + v ), c = f = i = —(w + w) 
Thus with v = a + if3,w = , -f + i5, the dephased matrix A is equal to: 
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Let us introduce now the following variables: 

x = —3a + (iw 6 — iw 3 )(3, y = —3a + (iw 3 — iw 6 )/3 
z = -37 + (iw 6 - iw 3 )5, t = -37 + (iw 3 - iw 6 )5 
Observe that we have in fact iw 6 — iw 3 = In terms of x, y, z, t, we have: 
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Thus we have reached again to the conclusion in the statement. 
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The general case similar, because when assuming A to be dephased, we obtain a p x p 
block decomposition of A as in the above examples, and this gives the result. For a full 
proof here, we refer to the proof of Proposition 4.5 below, at a = 2. □ 

Let us add now the trivial deformation part, and reformulate the result a bit as in 
Proposition 4.2 above. The statement that we obtain is as follows: 

Proposition 4.4. For F = F p 2, the elements A e T F C N are the solutions of 

A.. _ rOO , rOl , rlO , r 02 , r 20 , r 11 
A tJ — -^00 ~+~ ~+~ ~+~ ~+~ ^iO ~+~ Hj 

where the L variables are free, and dephased, and i,j are the reminders of i,j modulo p. 

Proof. By adding trivial deformations to the formula in Proposition 4.3, and then by 
reparametrizing these trivial deformations as in Proposition 4.2, we get: 

A^ = a + <2j + bj + M-fj 

Here Oj, bj are free variables, with the conventions ai = for % ^ 0, and bj = for j ^ 0. 
Now by using the further splitting of indices modulo p, we can write: 

Aij = a + fa + n + Si + Ej + M-q 

Here j3, 7, 5, e are now dephased, and this gives the formula in the statement. □ 

We are now in position of stating and proving our main technical result. We say that 
a matrix V s over the group Z p r x Z pS is dephased if its nonzero entries belong to: 

Here, and in what follows, we use the convention Z p -i = 0. 
Proposition 4.5. For F = F pa , the elements A e T F C^ are the solutions of 

A — \ T rs 

r+s<a 

where the L variables are free, and form dephased matrices L rs . 

Proof. Observe first that at a = 1, 2 this follows from Proposition 4.2 and Proposition 4.4 
respectively. Observe also that in the general case, the number of L variables is: 

d = ^ ^ I Zpr — Zpr — 1 J " J ZpS — ZpS — 1 J = ^ ^ p \^j p r — Z p r— 1| 

r+s<a r<a 
a 

= p a + J2p a ~ T '(P r - =P a + a (P - V 1 = (P + ap~ a)p a ^ 

Thus the number of L variables equals the defect d(F), so it is indeed the good one. 
As for the proof now, in the general case, this is quite similar to the one at a = 1, 2. 
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More precisely, consider the map L — > A. This map is linear, and in view of the above 
calculation, it is enough to prove that this map is injective, and has the correct target. 
For the injectivity part, recall that at a = 2 the formula in the statement reads: 

4 _ rOO , rOl , rlO , r02, r20, rll 
SHj — -^00 ^0,pj ^pifi "r ^iO ~*~ n pi,pj 

Now assume A = 0. Then with % = j = we get Lqq = 0. Using this, with % — 
and pj = 0, j 7^ we get Lqq + Lq| = 0, and so = 0. So, with i = and ^ we 
therefore obtain Z^gg + Lg| + L^ pj = 0, and so — 0. Now the same method gives as 
well succesively = and L™ = 0, so we are left with A^ = Lj^ pj , so we must have 
Lp}, P j = as well, and we are done. This method works of course for any a G N. 

Regarding now the "target" part, we must prove A G T f Cn- The equations are: 



E,„(H> ( V" 1 jrs _ jrs \ 

w I / j ^p a - r i,p a - s k ^p a - r j,p a - s k j 

k \r+s<a J 



= 



So, for any indices i,j and any r + s < a, we must prove that we have: 

E „,(»-.?> ( r rs — t ts \ — n 

V 1J p a - r i,p a ~ s k lj p a - r j,p a - s k) ~ u 

fe 

In order to do this, consider the following quantity: 



P k 

We must prove X^-j = Xj^-j. But, with k = m + p s n, we have: 

Y L \ ^ ln K m +P Sn ) j rs 

mn 

— \^ in lpSn \^ in lm T rs 
n m 

= SlO ^2 wlmL p*-ri,pa-s m 

m 

Thus we have I ^ =>• Xj/ = 0, and so = and we are done. □ 

Let us discuss now the affine deformations: 
Proposition 4.6. If N = p a , the defect of Fn is exhausted by affine deformations. 
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Proof. With k = m + p s n, as in the proof of Proposition 4.5 above, we have: 

k k r+s<a 

L rs _ L r 



n m r+s<a 



5 m p u JjP s m 

rrt r+s<a 

Now since this quantity vanishes for « 7^ j, this gives the result. □ 
Finally, let us collect all the above results into a single one: 
Theorem 4.7. For an isotypic Fourier matrix, H = F N with N = p a , we have 



lj ^ ^ p a-r itp a- Sj 



T° H C N = T H C N = T H C N ={Ae M N (R) A l3 = £ L rs 

r+s<a 

where the L variables are free, and form dephased matrices V s . 

Proof. This follows indeed from Proposition 4.5 and Proposition 4.6. □ 

Observe that this result shows that Conjecture 2.8 holds for the isotypic Fourier ma- 
trices. We will see in the next section that the same happens for any Fourier matrix. 

5. The general case 

In this section we discuss the computation of the enveloping tangent space T H C^ in 
the general case, where H = Fn is an arbitrary Fourier matrix. We will need: 

Lemma 5.1. If G = H x K is such that (\H\, \K\) = 1, the canonical inclusion 

Tf h C\h\ <S> T Fk C\k\ C T Fg C\g\ 
constructed in Proposition 3.1 above is an isomorphism. 

Proof. We have Fq = Fh x k, and the defect of this matrix is given by: 

x - \H x K\ x - I if x K\ 
d(F HxK )= ^ ^(h^) = 2. ord (h)ord(k) = 

{h,k)&HxK v ' ' (h,k)eHxK v ' v ' 

Thus the inclusion in the statement must be indeed an isomorphism. □ 

With this lemma in hand, the idea now will be simply to "glue" the various isotypic 
formulae coming from Theorem 4.7. Indeed, let us recall from there that in the isotypic 
case, N = p a , the parameter set for the enveloping tangent space is: 

X(p a ) = I I (Z p r - Zpr-l) X (Z p s - Z p s-l) 
r+s<a 
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Now since the defect is multiplicative over isotypic components, the parameter set in 
the general case, N = p" 1 . . .p a k k , will be simply given by: 

X(p a 1 \..pl*)=X(p?)x...xX(p a k ») 

One can obtain from this an even simpler description of the parameter set, just by 
expanding the product, and gluing the group components. Indeed: 

Definition 5.2. Given a finite abelian group G = Z p n x . . . x Z p ^ fc we set: 

G° = (Z p? - Z p? -0 x ... x (Z p? - Z p? -i) 

A matrix L e Mgxh(M) will be called dephased if L^- = for any G° x H° . 

Observe now that, with the above notation G°, the parameter set discussed above is 
given by the following simple formula: 

X(N) = |J G° x H° 

GxHCL N 

In addition, we can see that the collection of dephased matrices L e Mgxh(^) , over 
all possible configurations G x H C Z^, takes its parameters precisely in X(N). 
In order to formulate our main result, we will need one more definition: 

Definition 5.3. Given N = p^ 1 . . .p a k k and a subgroup G C Z N , we set 

<PG(ii,---,ik) = (Pi 1 ~ ri H,---pt k ~ rk ik) 
where the exponents r\ < ai are given by G = Z n x . . . x Z r k . 

Observe that in the case k — 1 this function is precisely the one appearing in Theorem 
4.7 above. In fact, we have the following generalization of Theorem 4.7: 



Theorem 5.4. For H = F N the vectors A e T H C^ appear as sums of type 

GH 

J <PG{i)<pH(j) 



a..- \" r GH 



where the L variables are free, and form dephased matrices L GH e Mqxh 

Proof. According to the above discussion, we just have to glue the various isotypic for- 
mulae coming from Theorem 4.7, by using Lemma 5.1. The gluing formula reads: 

A- • • • A- ■ A 



Ernsipi \ r t 

ai— ri . ai— si . • • • / ^ 



TkSkPk 



Pi 'lift' Jl ^ ' Pfc ^fc.Pfc fc Jfc 

vri+si<ai r k +s k <a k 



• • • / ai — ri . ai — si . • • • -Lj 



rkSkPk 



Pi' 'nft' 'Jl p fc fc fc «fe,P fe fc fc Jfc 
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Now, let us introduce the following variables: 



Tr 1 ...r k ,s 1 ...s k _ rnsi Tjk s k 
il...ik,jl—jk hjl ■ ■ ■ ikjk 



In terms of these new variables, the gluing formula reads: 

A ST^ rr 1 ...r k ,s 1 ...s k 

r 1 +s 1 <a 1 r k +s k <a k 

Now since the new L variables are obtained from the old L variables via a tensor product 
gluing, they are independent as well. Together with the fact that these variables form 
dephased matrices, in the sense of Definition 5.2 above, this gives the result. □ 

Corollary 5.5. The rationality conjecture holds for the Fourier matrices. 

Proof. Indeed, the formula in Theorem 5.4 shows that for H = the rational defect, as 
constructed in Definition 2.7, counts the same variables as the usual defect. □ 

A few comments now regarding the generalized Fourier matrix case, F = Fq. In the 
isotypic case G = TL^ x . . . x Z p a fc , according to Proposition 2.3 (2), we have: 



Wo) = £-L = £ 



ord(g) ^ max(or%i),...,or% fc )) 

ai <n ^ai+...+a fc 

I pmax(ri,...,r fc ) 



— • • • \%p r l — Zp^-i | . . . \ Z p r k — Z p r fc - 

ri=0 r k =0 

The combinatorics here is obviously much more complicated. In fact, if we assume 
a i < a 2 < • • • < a k then the defect formula from PQ , that we believe optimal, is: 

d(F G ) = N ( 1 + ^ p (*-r)«r-i+(«i+...+«r-i)-l( p *-H-l _ ^ _ — l]pfe- 
\ r=l 

Here ao = 0, and we use the standard notation [a] q = l + q + q 2 + ... + q a ~ x . 
As a conclusion, the following problem is open: 

Problem 5.6. How does TpC^ decompose, for a generalized Fourier matrix F = Fq? 

Regarding now the affine tangent cone TpCjq, we believe that, at least in the G = Z N 
case, an iterated application of Theorem 3.6 should give the answer, but we don't have a 
proof for this fact. As for the tangent cone TpC^, very little is known here, see [3]. 

Finally, in the general Butson matrix case, one technical problem comes from the reg- 
ularity conjecture in [2], mentioned in section 1 above. Indeed, a finer version of this 
conjecture would be that the regular complex Hadamard matrices are exactly the affine 
deformations of the Butson matrices, and this is of course very related to the present 
considerations. We intend to come back to these questions in some future work. 
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6. Probabilistic aspects 

We have seen in the previous sections that for certain complex Hadamard matrices the 
defect is the number of 1 entries, and the geometry basically comes from this. 

In general, the situation is of course much more complicated than that. One problem 
with the formula d(H) — |1 e H\ comes from the fact that the defect is insensitive to the 
equivalence relation in Definition 1.3, while the number of 1 entries is highly sensitive to 
it. So, unless H is given to us in some natural, standard form, as is the case for instance 
with the Fourier matrices, we have to take into account all the quantities of type |1 e H\, 
with H ranging over matrices which are equivalent to H . 

With this observation in mind, here is a first conjecture: 

Conjecture 6.1. We have the estimate d(H) < max |1 e H\, with the max ranging over 
all matrices H which are equivalent to H . 

Let us try now to formulate a finer version of this conjecture, roughly stating that 
u d(H) can be recaptured from the statistics of |1 e H\, over the matrices H ~ if". 

There is an obvious problem with this latter statement, coming from the fact that 
the number |1 e H\ is generically equal to 0. In order to overcome this issue, one idea 
is to restrict attention to the Butson matrices, and to allow in Definition 1.3 only the 
multiplication on rows and columns by the corresponding roots of unity. 

More precisely, let us denote Z s the group of s-th roots of unity, and by Cjy(s) the set 
of N x N complex Hadamard matrices having entries in Z s . We have then: 

Conjecture 6.2. For H e Cn(s), with s£N chosen to be minimal, we have 

min |1 G H\ < d(H) < max |1 e H\ 

with the min/max ranging over matrices of type = OibjHij, with a iy bj G Z s . 

Observe that in this statement we have dropped the action of the symmetric group on 
the rows and columns of H, because this action leaves invariant both the defect and the 
number of 1 entries. As for the assumption that s6N has to be minimal, this is of course 
in order for the lower bound to be non-trivial, because at s » this minimum is 0. 

Once again, it is not clear on how to approach this kind of statement, so instead of 
doing so, we will rather try now to further build on it, with some finer speculations and 
conjectures. So, let us go back to the comment following Conjecture 6.1. The quite vague 
statement formulated there can be now given a precise meaning, by using: 

Definition 6.3. Let H e Cn(s) be a Butson matrix. 

(1) We define if : Zf x Zf ->■ N by <p(a,b) = = 1}. 

(2) We let ji be the probability measure on N given by //({&}) = P((p = k). 
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In this definition P denotes the probability with respect to the uniform measure on the 
group Z^ x Z^ . In other words, we regard if as a random variable over this group, and 
we denote by \i the distribution of this random variable: 

1 



-,2N 



<f(a, b) = A; j 



#{(a,6)GZf xZf 

As a first observation, Conjecture 6.2 above can be reformulated as follows: 
Conjecture 6.4. For H G Cj^(s), with seN chosen to be minimal, we have 

d(H) G supp(fi) 

where the measure fi is the one constructed in Definition 6.3 above. 

Indeed, the fact that this statement is equivalent to Conjecture 6.2 follows from the 
fact that the image of the map ip constructed in Definition 6.3 consists exactly of the 
quantities |1 G H\ used in the formulation of Conjecture 6.2. 

Summarizing, we have reached to a quite conceptual reformulation and generalization 
of our very first statement, Conjecture 6.1 above, at least in the Butson matrix case. 

We will be back in a moment to this support problematics. But, let us formulate now 
yet another statement, which is our main conjecture on the subject: 

Conjecture 6.5 (Main conjecture). For H G Cn{s), with s£N chosen to be minimal, 
d(H) can be recaptured from the knowledge of the associated measure fi. 

As a first observation, this doesn't exactly generalize Conjecture 6.4. However, it is 
hard to imagine that a proof of this conjecture won't solve as well Conjecture 6.4. 

In order to further comment on this conjecture, let us first do some computations. The 
very first problem concerns of course the support of /i, and we have here: 

Proposition 6.6. The support of fi, with s G N chosen to be minimal, is as follows: 

(1) ForF 2 we get {1,3}. 

(2) For F 3 we get {0, 1,2,3,4,5}. 

(3) For F A we get {0, 1,2,3,4,5,6,7,8}. 

(4) For F 2j2 we get {4, 6, 8, 10, 12}. 

(5) ForF 5 we get {0,1,2,3,4,5,6,7,8,9}. 

Proof. Here is the proof for F4, in logarithmic form, with the operations being those in 
Definition 1.3, and with the subscripts denoting the total number of entries: 
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The proof for the other matrices in the statement is similar. □ 

Observe that Proposition 6.6 supports of course Conjecture 6.4, and that, in addition, 
it tells us many other things. Probably the most obvious one is that, in all the cases 
under consideration, the support is an arithmetic progression. We do not know how to 
prove this statement, so we would like to formulate this as a problem: 

Problem 6.7. Is it true that, for H G Cjy(s), with s G N chosen to be minimal, the 
support of n is always an arithmetic progression? 

Our first comment here is that the assumption that s is minimal cannot be dropped. 
Indeed, if we consider the matrix F 2i2 as an element of C 4 (4), and construct the measure 
fj, as in Definition 6.3, by using this exponent s = 4, the support that we obtain is: 

supp{fi) = {0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12} 

This can be indeed proved by proceeding as in the proof of Proposition 6.6. Now, 
observe that 11 is missing! So, we really need to assume that s is minimal. 
Let us try now to understand what happens at s = 2. Here we have: 

Proposition 6.8. Let H G Cn{2) with N > 2, and S C N be the support of fx. 

(1) S C 2N. 

(2) 5c[?,f]. 

(3) S is symmetric with respect to N 2 /2. 

Proof. Since we have iV > 2, we know that 4\N. 

(1) Let rj, Cj be the number of 1 entries on the rows and columns of H, so that we have 
|1 G H\ = Yl r i = c i- From the orthogonality condition for the rows and columns we 
see that the numbers rj have the same parity, and that the numbers q have as well the 
same partity. Indeed, let a,b,c,d G N be the lengts of the blocks in two distinct rows: 

1 1 -1 -1\ 
1-1 1 -I J 

We have then a + d = b + c = N/2 and a + b = r^a + c = rj, and since N/2 is even, 
we conclude that modulo 2 we have r j — r j = b — c = b + c = 0, as claimed. 

Now since is a sum of iV numbers having the same partity, the result follows. 
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(2) We pick two arbitrary rows of our matrix, say rows i,j, and use the study in (1). 
Since we have q = a + b and Cj = a + d we obtain: 



Ci + Cj b + c N 



a H = a H G 



N 3N 

T'T" 



4 

Now by averaging over all possible values of i ^ j, this gives the result. 

(3) This follows by switching signs over all the matrix. □ 

Observe that Proposition 6.8 is in fact very far from Problem 6.7. In fact, it is not even 
clear how to solve Problem 6.7 in the case of the Walsh matrices, and we have: 

Problem 6.9. What is supp(p) for a Walsh matrix, W N with N = 2 k ? 

According to the above results, at N = 2 we have supp(p) = {1, 3}, and at N = 4 we 
have supp(p) = {4,6,8,10,12}. It is also quite elementary to prove that at N = 8 we 
have supp(p) = {22, 24, ... , 40, 42}. In general, we do not know the answer. 

Let us discuss now the computation of the measure p itself. This appears to be a quite 
difficult task, and we have here only two results. The first one concerns the Fourier matrix 
F 2 . Here it is easy to see that p is uniformly distributed on {1,3}: 

\i= ^i + $0 

However, in view of the all the above issues with the minimality of s, it is perhaps 
interesting to see what happens for F 2 , at other values of s. And our result here is: 

Proposition 6.10. For F 2 G C 2 {s) with s even we have 

H = 4p* 3 - 6p* 2 + 4p-5 

where p = ^o~o + ^i is the rescaled spectral measure of the main character of Z s . 

Proof. We use the logarithmic writing. Consider the following matrix: 

% + a i + b 
J + a 3 + b + 1, 

Here the numbers i,j,a,b range in the set {0,1,..., s — 1}, and are taken modulo 
s. What we have to do is to examine the number of entries of F 2 , and compute the 
corresponding probability distribution p, which is supported on {0, 1,2,3}. 

A straightforward computation here gives the following formula: 

p = ^((s 3 - 4s 2 + 6s - 4)<5 + (4s 2 - 12s + 12)^ + (6s - 12)5 2 + 45 3 ) 

But this is the formula in the statement, and we are done. □ 

Summarizing, we have a complete picture of what happens for the first Walsh matrix, 
W 2 = F 2 . Our second result concerns the second Walsh matrix, W4 = F 2>2 . Here the 
computation at s G 2N arbitrary looks quite complicated, but at s = 2 we have: 
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Proposition 6.11. For the second Walsh matrix, i^,2 G C^(2), we have: 

fx = k^ + 1266 + 6(58 + 12Si ° + 5i2) 

Proof. It is technically convenient to use the equivalence — where the matrix .fQ 
is the one having —1 on the diagonal and 1 elsewhere. Let us perturb now this latter 
matrix, as in Definition 6.3. In logarithmic notation, we have: 

{Ki)ij = at + bj + Sij 

Thus if we want to compute the number of 1 entries, we have: 

\1 e K 4 \ = j,Oi = !>]} + #{(i,j)\i = j, at ^ty} 

= ^3, a i = b j} + ^ 

= #{(«, j)\a>i = bj} - #{i\ai = k} + #{i|a; ^ h} 

= #{(i,j)\a i = b j } + 4-2#{i\a i = b i } 

Now by writing the 16 x 16 tables for the two quantities appearing on the right, then 
making the above linear operation above, we obtain the explicit 16 x 16 table of the values 
of tp, and counting these values leads to the measure in the statement. □ 

Observe that for H e CV(2), computing the upper edge of the support of \i is the same 
as solving the corresponding Gale-Berlekamp game [H], [33]. So, we have: 

Problem 6.12 (Gale-Berlekamp game). Given a Butson matrix H e Cn{s), consider 
the matrices H obtained from it by multiplying the rows and columns by roots of unity of 
order s. What is the maximal number of 1 entries, over all these matrices H? 

In fact, the questions raised in the present section are related as well to the various 
algebraic or analytic problems regarding the 0-1 matrices, cf. |llj . [29] . [39J. 

Regarding now the fundamental quantum algebra/algebraic geometry problem men- 
tioned in the introduction, the previous considerations on |lj, involving quantum groups 
and Diaconis-Shahshahani type variables |12j . are still waiting to be explored. 

However, in view of the results found here, we believe that a better problem, at least for 
the moment, would be that consisting in connecting the above-mentioned Gale-Berlekamp 
type game to Jones' spin model problematics [19], [20], [21] . Indeed, the spin model 
problem is very related to the count of the 1 entries of the tensor powers H® k , and so can 
be regarded as "higher level" version of the Gale-Berlekamp type game. 

We intend to further explore this point of view in some future work. 
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